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ABSTRACT 

The stability of the two-dimensional flow induced by the tranverse 
oscillation of a cylinder in a viscous fluid is investigated in both the 
linear and weakly nonlinear regime. The major assumption that is made to 
simplify the problem is that the oscillation frequency is large in which case 
an unsteady boundary layer is set up on the cylinder. The basic flow induced 
by the motion of the cylinder depends on two spatial variables and is periodic 
in time. The stability analysis of this flow to axially periodic disturbances 
therefore leads to a partial differential system dependent on three variables. 
In the high frequency limit the linear stability problem can be reduced to a 
system dependent only on a radial variable and time. Furthermore, the 
coefficients of the differential operators in this system are periodic in time 
so that Floquet theory can be used to further reduce this system to a coupled 
infinite system of ordinary differential equations together with uncoupled 
homogeneous boundary conditions. The eigenvalues of this system are found 
numerically and predict instability entirely consistent with the experiments 
with circular cylinders performed by Honji [1981]. Results are given for 
cylinders of elliptic cross section and it is found that for any given 
eccentricity the most dangerous configuration is when the cylinder oscillates 
parallel to its minor axis. Some discussion of nonlinear effects is also 
given and for the circular cylinder it is shown that the steady streaming 
boundary layer of the basic flow is significantly altered by the instability. 
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INTRODUCTION 


Our concern is with the stability of a class of flows which exhibit the 
phenomenon usually referred to as "steady streaming." In particular, 
motivated by the recent experiments of Honji [1981], we consider in detail the 
stability of the flow induced by the transverse oscillations of a circular 
cylinder of radius a in a viscous fluid of kinematic viscosity v. This 
flow has been investigated by several authors following the boundary layer 
approach used by Schlichting [1932] . For a detailed discussion of the steady 
streaming induced by the oscillation of the cylinder the reader is referred to 
the papers of Stuart [1966] and Riley [1967]. 

The experiments of Honji [1981] illustrated clearly a phenomenon 
surprisingly not reported in previous experimental investigations of the flow. 
We refer to the observation made by Honji that the two-dimensional flow 
induced by the motion of the cylinder is unstable to axially-periodlc vortices 
of the Taylor-Gortler type at sufficiently large values of the amplitude of 
oscillation of the cylinder. The instability occurs in the Stokes layers at 
the cylinder in the locations where they are parallel to the direction of 
motion of the cylinder. The instability is apparently of the centrifugal type 
and is initially in the form of vortices aligned with the local flow 
direction. However, the steady streaming associated with the basic flow 
convects the dye used to visualize the vortices away from the Stokes layer. 
At larger amplitudes of osciallatlon the dye streakes produced by the vortices 
disappear and the flow is said by Honji to be turbulent and separated. It was 
suggested by Honji that the instability might be of the type which is known to 
occur in a Stokes layer on a torslonally oscillating cylinder and it is this 
possibility that we shall investigate in this paper. 
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This latter type of instability has been investigated in detail by 
Semlnara and Hall [1976, 1977], Park, Barenghi and Donnelly [1980] and Hall 
[1981], Suppose then that an infinitely long cylinder oscillates torslonally 
about its axis with angular velocity Q cos wt in a viscous fluid. At 
sufficiently small values of w the flow is purely circumferential whilst if 
0 ) is slowly Increased then at a critical value of w = an array of 

vortices, periodic along the cylinders, develops in the boundary layer at the 
cylinder. The strength of the vortices increases as w is increased further 
but at a second critical value of w, say '^^ 2 * vortices Interact with 

each other and the flow rapidly becomes turbulent with no apparent periodicity 
along the cylinders. The theoretical description of this flow for w < 
given by Semlnara and Hall [1976, 1977] was verified experimentally by Park, 
Barenghi, and Donnelly [1980]. However, the secondary stage of the 
instability is perhaps only partially explained by the subharmonic instability 
mechanism described by Hall [1981]. It is of Interest to note that the values 
of ‘^(.2 quite close so that, once the initial vortex structure has 

been set up, a relatively small Increase in m leads to a turbulent flow. 
Thus, if this instability mechanism is Indeed operating in the experiments of 
Honji, it is clear that the two-dimensional flows of the type discussed by 
Schllchting [1932], Stuart [1966], etc. will be greatly altered. The primary 
aim of the present investigation is to determine the parameter range in which 
the two-dimensional solution is a stable solution of the Navler-Stokes 
equations. The feature of the basic flow which makes a stability calculation 
nontrivial is, of course, the fact that the basic flow depends upon time and 
two spatial coordinates. 

It is to be expected that at sufficiently large amplitudes of oscillation 
the boundary layer at the cylinder will separate and the attached flow 
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strategy of the type discussed by Schlichting and subsequent authors fails. 
The results of Honji suggest that this does not occur before the instability 
mechanism is operational. Thus our stability calculation provides an upper 
limit for the oscillation amplitude beyond which there is no reason to compute 
the basic flow. It is therefore likely that any laminar separation theory for 
the oscillating cylinder problem will not be relevant to experimental 
observations. However, the success of laminar separation theory in high 
Reynolds number steady flows (see Smith [1982]) suggests that this might not 
be the case. 

Suppose then that a circular cylinder of radius a oscillates with 
velocity Uq cos wt along a diameter in a fluid of viscosity v. The 
parameters which govern the two-dimensional flow are 



(1.1a) 


X 


toa 


R 

s 



(1.1b) 

(1.1c) 


The frequency parameter P is taken to be large so that the unsteady boundary 
layer on the cylinder is thin compared to its radius. The parameter X 

represents the ratio of the amplitude of oscillation of the cylinder to the 
cylinder radius and is taken to be small. Stuart [1966] has discussed the 
crucial role played by the steady streaming Reynolds number Rg in 

determining the nature of the steady streaming set up outside the Stokes layer 
on the cylinder. If Rg is small the motion is determined by solving the 
Stokes equations whereas for large Rg an outer boundary layer flow exists. 



-4- 


In order to obtain some idea about the parameter regime in which an 

instability might occur we note that the first order oscillatory flow set up 

by the motion of the cylinder is confined to a thin layer of thickness 
, I/9 

(v/m) ^ at the cylinder and so the radius of curvature of the paths of fluid 

particles is of order a. Thus, the Taylor number which characterizes this 


u: 


boundary layer flow is of order 


0 


1/2 3/2 ^s ^ 


I/2 


The instability 


a V u) 

mechanism described by Seminara and Hall [1976, 1977] operates when this 
Taylor number is 0(1) so we conclude that in the present problem the regime 

A V2 


of interest is R 


For this reason we confine our attention in this 


s 

paper to the stability of the two-dimensional flow around the cylinder in the 

V 9 

limit 6 -> 0 °, R ~ 0(3 ^ • We further note that in this limit A is 
s 

- 1 /, 

0(3 so that the boundary layer on the cylinder is essentially a Stokes 


layer. 

The above comparison between the torsionally and transversely oscillating 
cylinder flows Ignores the spatial variation around the cylinder of the first 
order boundary layer flow in the latter case. An examination of this 
structure shows that the flow is locally most unstable at the positions 
6 = ± ir/2 if the direction of oscillation is along the x-axis. We shall show 
that a self-consistent as 3 miptotlc description of the linear stability problem 

I/9 

is possible for R ~ 3 , 3 Furthermore, we show that the Instability 

“1 /8 

is confined to 3 neighbourhoods of the positions 6 = ± ir/2. More 

precisely we show that the flow is formally unstable when 


R > R 
s sc 


= Rn 5 


V2 


+ R, 



+ Ro e 


1/8 




( 1 . 2 ) 


where Rq , Rj^ , 


etc. are 0(1) constants to be evaluated. 


In fact, we 
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determlne only Rq, Rj^ and find that the resulting critical value of X 

agrees almost exactly with the experimental results of Honji. 

The analysis used for the circular cylinder problem can be easily modified 
to more complicated steady streaming flows. We shall show that, having 
investigated the linear stability problem for the circular cylinder we can, to 
first order in 3 , write down the critical stability parameter with only a 
knowledge of the first order outer potential flow. However, at the next order 
there are technical differences between the circular cylinder problem and for 
example the problem associated with elliptic cylinders. We shall see that 
these technical differences depend on whether or not the stagnation point of 
attachment of the steady streaming coincides with the most unstable part of 
the boundary layer. 

Va 

Some discussion of nonlinear effects for {R - R } ~ 0(3 is given. 

s sc 

For the circular cylinder problem a strong interaction between the steady 

streaming and the instability occurs. In fact it appears that the higher 

modes of instability lead to the separation of the steady streaming boundary 

—1 /8 

layer within an angle 0(3 ' ) of the point of attachment of the layer. 

The procedure adopted in the rest of this paper is as follows. In Section 

lA 

2 the linear stability problems is formulated for 3 ^ R^ ~ 3 ^ and an 
as 3 unptotic solution of the problem is given. In Section 3 the results of the 
numerical solution of the eigenvalue problem obtained in Section 2 are given 
and compared to Honji 's experimental observation. 

In Section 4 we discuss the relevance of our calculations to more 
complicated flows. More precisely we consider the stability of the flow 
induced by the oscillation of an elliptic cylinder. We consider an ellipse 
with major and minor axis a,b with major axis inclined at angle a to the 
x-axls in which direction the cylinder is oscillating. We find that. 



depending on the values of a and b/a, there are two or six locations where 
instability will occur. In Section 5 we consider the nonlinear development of 
the instability whilst in Section 6 we draw some conclusions. 


2. FORMULATION AND SOLUTION OF THE LINEAR STABILITY PROBLEM IN THE LIMIT 3 + » 
The first step in our formulation is to note that by a simple change of 
axes we can take the cylinder to be held fixed whilst the fluid at Infinity 
oscillates with speed Uq cos o)t parallel to the x-axis. It is convenient 
for us to work in cylindrical polar coordinates (r,6,z') with the z' axis 
along the axis of the cylinder. We now define the variables n , z, and t by 


n = (r-a)(^) ’■^2 ^ (2.1a) 

z = (2.1b) 


T = 0)t. 


(2.1c) 


Following the scalings discussed in Section 1 we write 


R 

s 


_ T3 ^^2 

p/2 


( 2 . 2 ) 


where T is 0(3*^) and is, of course, the Taylor number. We shall 
investigate the stability of the boundary layer on the cylinder in which the 
basic velocity field is (u,v,0) with 
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V = 2Uq{s1ii 0 VQ(n,x) + v^(ti,t) + ••• } 


(2.3b) 


where in particular 


Vq = cos T - cos(x-p)e ^ , 


whilst Uq, etc. can be found in, for example, Stuart [1966]. 

Following the scalings used by Seminara and Hall [1976] we perturb the 
basic flow such that the new velocity field is 

u = (u,v,0) + (/ 2vu) U(n,0,z,x), UpV(ri ,0 ,z ,T ) , / 2vu W(ri ,0 ,z,x)) , (2.4) 

whilst the corresponding pressure perturbation is puv P(ri,0,z,x). If the 
above expression is substituted into the Navler-Stokes equations we find 
that U, V, W, and P satisfy 


L'U = - 2T sin 0 Vq V + + 0(0 ), 


(2.5a) 


L'v = T'^ T 


3/8 V 2 p 3/4 |P ^ Q 0 y ^ 1/4)^ (2.5b) 


L'W = ||- + Q3 + 0(0 ), 

H 4 . 2~^^^ ^ ^ ^ 

3n p V 4 90 9z 


(2.5c) 

(2.5d) 


Here the nonlinear terms Q;]^, Q 2 , and Q 3 are given by 
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whilst the operator L' has been defined by 

8z^ 


( 2 . 6 ) 


We further note that the 0(3 ) terms not shown explicitly in (2.5) 
comprise both linear and nonlinear terms. However the linear terms not shown 
explicitly vanish when 0 = ± tt/2 and for that reason are negligible in the 
following analysis. The nonlinear terms not shown explicitly do not vanish at 
0 = ± tt/2 but the smallness of the disturbance which we assume in Section 5 
means that, to the order considered in this paper, these terms are also 
negligible. 

For the remainder of this section we neglect the nonlinear terms in (2.5) 
and assume that P, U, V are proportional to cos kz whilst W is 
proportional to sin kz. Here k is a constant axial wavenumber and it is 
now convenient to eliminate W, P from the linearized form of (2.5) to give 


- k^)u = 2k^ T sin 0 Vq V - 
9n 


2^^^ sin 0 v„ . _i, 

-tj~— 


9^0 - V/ 

LV = 4 sin 0 U + 0(3 , 


(2.7b) 
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which are to be solved subject to 


U = V = 


3U 

9ri 


= 0 , 


U,V > 0, 


n = 0 




( 2 . 8 ) 


The operator L appearing in (2.7) is simply L' with — j replaced by 

3z 

-k2. 

It can be seen from (2.7) that the 0 variation of the disturbance is 
slow compared to the x and n variations. The 0 dependence of U and 
V can therefore be taken care of by a WKB type of approach. However, since 
we are interested in the most unstable disturbances, it is convenient for us 
to use a multiple scale method. We can see from (2.7) that, ignoring the term 


proportional to 0 


V 4 0 


30 


, the "effective" Taylor number of the flow Is 


T sin" 0 which has local maxima at 0 = i: ir/2. Hence in the neighbourhood 
of say 0 = -rr/2 the effective Taylor number is t{1 - (0 - tt/2) + • • • } . 
The symmetry of (2.7) about 0 = ir/2 means that when the WKB formulation is 
used the point 0 = ir/2 is a turning point. Since the local Taylor number 
has a maximum at 0 = Tr/2 the turning point is of the second order and the 
usual scaling analysis shows that a transition layer of thickness 0(P ^^^) 

exists near 0 = tt/2. This situation is similar to that found by Hall [1982] 
who investigated the growth of small wavelength Gortler vortices in boundary 
layers on concave walls. In that problem the most unstable modes have a 
vertical structure concentrated in an internal transition layer which again 
corresponds to a second order turning point. 

The discussion above clearly applies to the neighborhood of 0 = - ti/2 
but let us concentrate on the transition layer at 0 = ir/2 and write 



10 


$ = [0 - tt/2] 


, 1/8 


We seek a solution of (2.7) by expanding U in the form 

U = Uq(p,t,$) + e~ Uj(n,T.$) + U2(n,x,$) + ••• , (2.8) 

together with a similar expansion for V. The Taylor number T then expands 
as 

T = Tq + T^ + 3" T^ + ••• . (2.9) 


Here we have anticipated that the first-order correction to T from Tq is 
0(3 ) rather than 0(3 ^^^) as might be expected from (2.8). The above 

expansions are then substituted into (2.7) and 0 replaced by 
tt/ 2 + 3 If terms of 0(3^) are equated we obtain 


•8 


8n 

„2 


k - 


- 


30 


^ 3 t 2 “ ^ ^^0 ~ '^0 ^0 ^0 

3p 


a 3v» 

2 — Iv = 4 - U 

3t 1'^0 3n ^0» 


(2.10) 


and the appropriate boundary conditions are 


8U 

^0 = ^ = ^0 = ^ 


Uq.Vq > 0, 


0 ^ 


( 2 . 11 ) 


The partial differential system (2.10), (2.11) governs the centrifugal 

instability of a Stokes layer on a cylinder driven by a pressure gradient 
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rather than by the motion of the cylinder as was the case in Seminars and Hall 
[1976],. The major difference is that in the latter paper the function Vq is 
replaced by cos(x-Ti)e . In order to determine the value of Tq = TQ(k) 
above which exponentially growing solutions of (2.10), (2.11) exist we seek 
periodic solutions of the latter system by writing 


Uq = A($) I Uj(n)e^'^^, Vq = A($) I v5(n)e^"^. 


The sequences of functions {Uq} , {v”} therefore satisfy the ordinary 
differential system 



dn^ 


u 


n 



d n 
'dp 0 


2in] [-^-y 
dn 

21n]v^ = 


- k^]u” = k^ TQ{[l-e"^^^"''^^]VQ"^ + [ 1-e”’^ 
-2{(l+l)e”’^^^'*"^^ + (l-l)e”’^^^”^^U^^} (2.12) 

p = 0 


Uq, Vq •> 0, n with n = 0,dbl,±2,»*» . 

The numerical solution of (2.12) will be discussed later and it suffices to 

say that the eigenrelatlon k = ^(Tq) can be determined. The amplitude 

function A($) remains undetermined at this order. 

—1 /8 

At order 3 the partial differential equations satisfied by 


are found to be 
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^ IrHfr - vj 

9n 


= 1^''^ A 8 [vJ^ - k^)u_ - V UJ 
0 9$ 0'' 2 ''0 onn QJ 

dri 


(2.13a) 


u2-2^W -«!!2o ^ !!» 

9tJ 1 9n 1 0 0 90 » 


together with boundary conditions identical to (2.11). The differential 
system for (Uj^, V^^) is an inhomogeneous form of (2.10), (2.11) and so we 
require that a solvability condition must be applied if the system is to have 
a solution. However, it can be inferred from (2.12) that the sequence of 
eigenfunctions ^^n^ such that either 


i) = 0, n even, ^n “ 


b) = 0, n odd, = 0, n 


In fact, our calculation showed that the most unstable mode corresponds to b) 
above. In either case we see that the Inhomogeneous terms on the right-hand 

il.T 

sides of the equations for and are proportional to e Uq, 

i i T 

e Vq in which case the solvability condition is automatically satisfied. 
The solution of the system for (U^, V^) can be written 




QA r inx „n . V 


(2.13b) 


inx „n 
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together with a similar expression for Here B is another amplitude 

function to be determined at higher order. 

_ 1 /, 

At order B ^ the function pair (U 2 ,V 2 ) is found to satisfy 


2 IrKfr - ’’o ''0 

8n 


3n 


V, 


= 2^/^ 1_ I- [9 


0 8$ ^ QL 2 

3n 


Uj) +2k^ v„(Tj_ 


T„ V, 


9ri 


- 2 


i-l V - 4!^0 - 2=/'' T^/2 V - 

9tJ 2 8n 2 0 9$ 


^^0 2 

- 2 - U 9 

9n ’ 


whilst the boundary conditions are again identical to (2.11). The forcing 
terms on the right-hand side of the above equations are synchronous with the 
solutions of the homogeneous forms of the equations and a solution for (U 2 ,V 2 ) 
will not in general exist. However, by considering the partial differential 
system adjoint to (2.10), (2.11), we find that a solution exists if 


d^A 

d$^ 


+ p[t^ - $^]a = 0 


(2.14) 


where p is given by 


2 r . + 


k / / u V„ v„ dn dT 


0 0 


'0 0 


V4.V2f J , 

0 0 


I,+ r, , — r,+ f 3 

U U, + v„ U — K- 
oqq 1 0 '-.2 

9p 


- dn dx 


(2.15) 


4* 4* 

and (U ,V ) satisfy the adjoint differential system 
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+ 2 IV-^ - 2k^ T„ 


;9 i2,-)9ir8 ,2i + , — + 

^TT" ^ +2 — }{ - k }U = 4 V V . 

9n 9n 


u+ = |IL = v+ = 0 . 

3ri ’ 


n = 0, 


u'^, v'^ -V 0, 


The eigenvalues of the adjoint system are of course identical to those of 
(2.10), (2.11) and can be obtained by Fourier expanding (U'*', V"^) to obtain 
an infinite set of coupled ordinary differential equations. 

The constant y is a function of k and our computation suggests that it 
is positive near the critical value of Tq. Such a result is not unexpected 
since it means that nonaxisymmetric disturbances to a Stokes layer on a 
cylinder are more stable than axis 3 mimetirc disturbances. This result was 
found by Duck and Hall [1980] for the case when the flow is driven by the 
motion of the cylinder whilst the present results show that this is also the 
case if a pressure gradient is driving the Stokes layer. The amplitude 
equation (2.14) has solutions which decay to zero when $ > if y > 0. 
These solutions are 


A($) = A^($) = U^(-n -1/2 ,2y 


(2.16) 


where is the nth parabolic cylinder function and the value of Tj^ 


corresponding to A^^ is 


T = T, =2 
1 In 


[n + 1/2 ] 


(2.17) 
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The least stable mode corresponds to 

Aq($) = e 


n = 0 in which case 
-(p ^^2 $2/2) 


(2.18) 


The expansion procedure described above can be continued to any order and we 
note here that the next nonzero term in the expansion of T is T2« 


3. m NUMERICAL SOLUTION OF THE LINEAR EIGENVALUE PROBLEM 

The solution of eigenvalue problems such as (2.12) is now a routine 
procedure and we shall give only the essential details of the calculations. 
The first step is to reduce (2.12) to a finite set of equations by setting 

= 0 for In I > M. We then replace «= by so that (2.12) has 

been approximated by a finite system of equations on a finite interval. Of 
course, it is necessary to vary and M to find appropriate values which 

enable us to solve (2.12) with sufficient accuracy. 

If p is sufficiently large then u”, v” satisfy 

rd^ ,2 ird^ ,2i„n ,2 „ r n-1 , ..n+li 

- k -2inJ[— ^-kJU =k Tq[Vq + Vq J, 

dp dp 0 

[-^ - k^ - 2in]VQ = 0. 
dp 

There are three independent solutions of these equations which decay 
exponentially to zero when p •> «>. Thus, if n is restricted to the range 
-M < n < M we can use these solutions to integrate the differential equations 
from h = P^j to p = 0 thus obtaining 6N + 3 independent solutions of the 
reduced system of equations. This integration was carried out using a fourth- 
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order Runge-Kutta scheme with step length h. These Independent solutions of 
(2.12) can be combined at n = 0 to satisfy 6M + 2 of the required boundary 
conditions there. The remaining boundary condition is automatically satisfied 
if k = k(TQ) is an eigenvalue of the reduced system. This eigenvalue of 
the reduced system will depend on M, and h but by Increasing and 

decreasing h an eigenvalue of (2.12) can be obtained. 

In our calculations it was found that M = 6, n = 10, and h =0.25 

gave results correct to the accuracy given in this section. In Figure 1 we 
have shown the neutral curve k = k(TQ) and the minimum of this curve 

corresponds to 

Tq = Tq^ = 11.99, (3.1a) 

k = k = .51 (3.1b) 

c 

It is interesting to note that the corresponding value of Tq for a 

torslonally oscillating cylinder in ~ 230. so that transverse oscillations of 
the cylinder produce a much more unstable flow. There is no obvious physical 
reason why this should be the case. 

The eigenfunctions corresponding to Figure 1 were normalized by taking 
Uq (0) =1 and have the property that ~ 0, = 0 when n is an odd 

integer. The functions U^, corresponding to the critical case are shown 

in Figure 2. We note that the disturbance is most pronounced near n ~ 3. 
It is interesting to note that an asymptotic solution of (2.12) in the limit 
a 00 with T ~ a^ shows that the vortices become concentrated in an 

internal viscous layer of thickness a" ^^2 . in this layer the functions 
Uq, Vq can all be expressed in terms of parabolic cylinder functions. A 
similar calculation for the torslonally oscillating cylinder problem shows 
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4 

that when a •> with T ~ a the vortices become concentrated in a layer of 
thickness a"^^^ near n = 0 and the eigenfunctions are then determined in 
terms of the Airy function Aj|^(x). 

In order to check the eigenvalues shown in Figure 1 the solution of the 
adjoint system was computed in a similar manner. The adjoint eigenfunctions 
were found to have the property U^, = 0 for n an odd Integer. The 

A A 

Inhomogenous system for was found by a shooting procedure similar to 
that used for (2.12). The integrals appearing in the definition of p were 
then evaluated using Simpson's rule. We obtained p = .033 with TQ,a as 
given by (3.1a,b). The critical value of Tj^ is therefore given by 

T, = 5.51 
Ic 

so that the critical value of Rg is 

R = R = 4.24[e^''2 + .46 + ... ]. (3.2) 
s sc 

If Rg is greater than Rg^ the vortices grow exponentially in time but 
remain localized near 6 = ir/2. In order to compare our result with those of 
Honji we rewrite (3.2) in the form 

We further note that, in the notation of Honji, X is equal to the ratio of 
the cylinder oscillation amplitude dQ to the diameter D = 2a and that the 
Strouhal number St defined by Honji is related to 3 by 
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2 

St = - 8. 
n 

In Figure 3 we have compared our theoretical prediction of A = d^/D with 
Honjl's results. We recall that above the lower of the two sequences of 
experimental points Honji observed Taylor-Gortler vortices. There seems 
little doubt that the instability mechanism discussed here is responsible for 
the vortices seen by Hohji. Surprisingly we see that (3.3) is in excellent 
argument with Honjl's results even for X ~ 1 even though (3.3) is formally 
valid only in the limit 3 ^ 


4. LINEAR THEORY FOR MORE GENERAL STEADY STREAMING FLOWS 

We shall in this section discuss the modifications to the expansion 
procedure of Section 3 which are necessary when the basic flow does not have 
the symmetry of the circular cylinder problem. Suppose then that we consider 
the stability of the boundary layer induced by the outer potential flow 
Uq U(x) cos wt interacting with a rigid wall of local radius of (convex) 
curvature aR(x). Here x is a dimensionless variable which measures 
distance along the wall. We again take p to be a normal variable scaled on 
the Stokes layer length scale (Zv/oi) ^ and the basic flow in the boundary 
layer will be of the form (2.3) with sin 0, cos 0, and sin 20 replaced by 
U(x), U'(x) and V 2 ^(x) U'(x) respectively. 

At any local station x along the wall the local Taylor number T^ which 
governs the stability of the boundary layer is given by 


= 


2^^^ Uq U^(x) 


a V 




R(x) 
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and instability of the localized type discussed in Section 3 will occur at 

2 -1 

X = which is a maximum of the function U (x) R (x). In the 

neighborhood of x^^^ we write 


$ = (x - X )B 
m 


+1/8 


and expand the disturbance as in (2.8). The only essential differences are 
that the coefficients of the amplitude equation (2.14) are altered because: 
(1) the radius of curvature must be expanded locally, and (2) the terms 

- V/ - Va 

corresponding to the 3 ^ terms in (2.3) remain 0(3 ^) near x = x^^^ and 

so contribute to the linear term in (2.14). However, the first-order term in 
the expansion of the Taylor number does not depend on these higher order 
alterations and so we can say that, correct to first order, the boundary layer 
is locally neutrally stable at x = x^^^ if 


d 

dx 




and T. = 11.99, x = x . 
i m 
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where Max(^) denotes the largest maximum value of ^ . We see that, 
having performed the calculation of Section 3, the critical configuration 
defined by (4.1) can to first order be written down with a knowledge of only 
the outer potential flow and the local radius of curvature. We Illustrate the 
simplicity of the procedure by considering the stability of the flow Induced 
by the oscillation of an elliptic cylinder in a viscous fluid. 

The basic flow induced by the high frequency oscillation of an elliptical 
cylinder has been given by Davidson and Riley [1972]. We suppose that the 

ellipse has major and minor axis of length 2a, 2b respectively and that the 

cylinder oscillates with velocity Uq cos o)t in a direction making an angle 
-a with the x axis. In order to find the critical Taylor number for the 

flow we require the first order potential flow and the local radius of 

curvature of the ellipse. The potential flow is obtained in a routine manner 
by mapping the ellipse onto a circle of radius V 2 (a+b) . If we use the 
parmetric representation of the ellipse 

X = a cos (|) , y = b sin o < <)) < 2ir , 

then the slip velocity of the potential flow is 

U.(l+K)sin((|)-a) 

U„U((|)) = ^ » where K = b/a. 

° (sin^ <[, + cos^ (f) 

The radius of curvature is easily obtained in terms of ()> and we find that 
the local Taylor number T^ is given by 

3/2 2 22 

2 ' Uq( 1+K) sin (())-a) 

V 2 3/2, , 2 ^ ^ ^2 2 ^,5/2 ’ 

av (D (sin ((> + K cos <j)) 
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and it follox^/s that the flow is locally neutrally stable at ~ 


d(j) 


0 , 


T„ = 11.99, 



We shall take a to he held fixed whilst K and a the angle of attack vary 
and compare the critical Taylor number T for the flow with that appropriate 
to a circular cylinder of radius a. 

The values of T^^ at which the flow is locally neutrally stable at some 
value of are found by considering the maxima of the function 


S((t)) = sin (({)-a)(sin <|) + K cos (j>) 


5/2 


(4.2) 


Unlike the circular cylinder case there can be more than two values of at 

m 

which this function has a maximum. In view of the symmetry of the problem we 

can restrict a to the range 0 < a < 7r/2 and (J) to the range 0 < <t> < ir . 

After some calculations it can be seen that the results for arbitrary 
values of a can be understood by first considering the case a = 0. In 
this case the function S(([)) has a maximum at (j> = (J)^ = ir/2 for all values 

of K. However, when K < /3/5 two further maxima occur at 


—1 I 2K^ 

(}) = <{)„ = sin / 7 f- and cj) = <{>„ = it - (j>„ . These two extra maxima 

3(1-K ) 

emanate from <)) == 0 when K = /3/5 and 

S((|)2) = S(<j>3) > S((j)j^), for K < /ITS . 

Hence for K < / 3/5 there are six potentially unstable points on the cylinder 
and the most unstable points are (j) = <\> 2 > <1)3, <1>2 » <)>3 + The migration 
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of the most unstable point away from 4 = 12 Is due to the relatively large 

increase in curvature away from ^ = tt/2 caused by increasing the 

eccentricity of the ellipse. We note that when a = 0, K = / 3/5 the scaling 

of Section 3 needs a more significant alteration since tj) = ir/2 is then a 

- 1/12 

fourth-order turning point. We must then choose to work in a 3 

neighbourhood of ir/2. The appropriate amplitude equation then has the term 

2 4 

proportional to $ A in (2.3) replaced by $ A. 

In Figure 4 we have shown the dependence of T/Tq (where Tq is the 

critical value for the case K = 1) on K with a = 0. We obtain a familiar 

cusp-shaped curve and we note that for K > /3/5, T/Tq is a single valued 

function of K. The lower curve for K < /3/5 corresponds to the two equally 

unstable locations <1)2 and (j)^ and passes through the origin. This means 

that the critical value of T/Tq can be made arbitrarily small by taking the 

limit K > 0. We can see in Figure 5 that in this limit the locations of the 

most unstable positions approach <t> = 0,ir where the radius of curvature is 

clearly greatest. 

The results for a 0 are obtained by describing the unfolding of 

Figures 4 and 5 when 0 < a « 1. The upper and lower curve to the left of 

the cusp move up and to the left when a increases from zero. The lower 

curve for a = 0 is in fact two coincident curves corresponding to ^ 

and (f) = The other one of these curves remains connected to (t/Tq = 1, 

K = 1) and (T/Tq, K = 0) but moves downwards until T/Tq Is eventually a 

monotonically increasing function of K on this branch. Ultimately the 

T 4K^ 

branch as3nnptotes to the curve j which corresponds to a = tt/2 

0 (1+K) 

whilst the detached upper branches rapidly move to be left and upwards when 
a increases. Finally when a = ir/2 there are only two maxima on the 
cylinder at <j) = 0 ,tt. In Figures 4 and 5 we have illustrated this process 
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for a few values of a. The curves I, II, and III of these two figures 
correspond. 

Suppose now that we have an elliptical cylinder with K fixed and we 

require the most stable or unstable orientation of this ellipse in an 

oscillatory flow. It follows from Figure 4 that if we wish to keep the flow 

stable then we choose a = 0 whilst if we wish to set up an unstable flow 

then we take a = tt/2. Next suppose that a the angle of attack is held 
fixed and K can be varied. We see from Figure 4 that for some values of a 
there is range of values for K which give a flow more stable than that 
around a circular cylinder of radius a. The most pronounced effect of 
Increasing eccentricity corresponds to the a = tt/2 case. Here we see that 
changing say K from 1 to V 2 produces a decrease in the critical Taylor 
number by a factor of ~20. 

5. THE NONLIMAR DEVELOPMENT OF THE INSTABILITY FOR THE CIRCULAR CYLINDER PROBLEM 


We shall now describe the manner in which finite amplitude effects become 
important close to the critical Taylor number. We suppose that Rg differs 
from its critical value by 0(0 ). On the basis of weakly nonlinear 
stability theory we expect that the azimuthal velocity component set up will 
be of magnitude 0(g^^^ The most illuminating choice for 6 is that 
which ensures that the azimuthal structure of the disturbance is determined at 
the same order as is the amplitude of the disturbance. The appropriate choice 
of 6 is 6 = V 4 and we therefore expand T in the form 


T = T» + 0 T + » • • 

UC 


where Tq^ is the critical value of Tq« 


(5.1) 
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The disturbance quantities U, V, W, and P then expand as 


U = U_ cos k z + [U, cos k + U- cos 2k zl 

0 c ‘•1 c 2 c 


rV4 


+ 


-3/8, 


U_ cos k z -f U, cos 2 k z -b cos 3k z + U ,_1 -I- 

■s r* A z' S r* MfiJ 


V = 3 V. cos k z -f 3 [v^ cos k z -I- V_ cos 2 k z + V.,«l 

0 c ^1 c 2 c MQJ 


+ 3 cos k z + V, cos 2k z -b cos 3k z -b V„, 1 + 

‘-1 lx n S n Ml J 


W = 3“’''^® W- sin k z + 3 [w, sin k z -b W_ sin 2 k zl 

0 c '-1 c 2 c ■’ 


+ 3 sin k z + W, sin 2k z + W,. sin 3 k z] -b 

3 c 4 c 5 c 


P = 3"^^^ P^ cos k z + [P, cos k z -b P. cos 2k z] 

0 c ‘•1 c 2 c ■' 


V4 


+ 3 ^^^fPo cos k z + P, cos 2k z -b P^ cos 3k z] 
3 c 4 c 5 c 


+ ■'mo + e'"'® ”«1 + 8 "'*'’m2 + 6-^'® P„3 + 


where apart from Pjyio* ^M1 > ^M2 which depend only on t and 

coefficients in these expansions are functions of t, $ and n . 

functions P^O’ ^Ml» ^M2 essentially pressure eigenfunctions 

to satisfy all the required conditions on the mean velocity field. 
DiPrima and Stuart [1975]) for a discussion on the need for 
eigenfunctions in centrifugal instability problems.) 


(5.2a) 


(5.2b) 


(5.2c) 


(5. 2d) 

$ , the 
The 
needed 
(See 
such 
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It is now a straightforward procedure to substitute from (5.2) into (2.5) 

— 1 /8 — 1 /8 

and successively equate like powers of 3 .At order 3 we find that 

(Uq, Vq) satisfy the linear stability problem (2.10) with Tq = Tq^, k = k^, 
so that 

(Uq.Vq) = A($)(Uq.Vq) , 



Here the functions {Uq} , {v^} satisfy (2.12) with k = k^,, Tq = Tq^. At 

- V/ 

order 3 ^ we find that the term proportional to sin k^ z, cos k^ z again 

satisfy (2.13a) so that (U^jVj^) is given by (2.13b). In addition to the 
fundamental modes generated at this order there are first harmonic and mean 
flow correction terms produced by the nonlinear terms Q 2 , and Qq 

appearing in (2.5). After some manipulation we can show that (U 2 ,V 2 ) and 
(UmQjVjjq) can be expressed in the form 


(Uj.Vj) - A^(«)(U2(n,T),V2(tl,T) 




A A 


where ( 02 ,^ 2 ) satisfy 
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r9^ ,, 2 „ 8 1 r9^ ..2i 


- 4k - 2 — y - 4k }U, - 8k T v_ V„ 

\ 2 c 9 t^ ~ 2 2 c Oc 0 2 

9ti 8ri 


/V A 


k^ T_ vj + 2(U- U_ - U- U. } , 
c Oc 0 '0 Oiiriil On Onn^ 


O ^ A A A A 

- 4k^ - 2 ^}V2 - 4 - V, = {U^ V^} . (5.3) 


9n 


^2 = ^2 = ^ ^2 = 0 . 


n = 0 


U2.V2 - 0, 


n ->• “. 


The mean flow correction terms Vjjq satisfies 


^2 9V„„ . 

1— (V ) - 2 — ^ = - 
. 2 '■ MO'' 9t 
9n 


8n ("0 'ol • 




(5.4) 


The forcing term on the right-hand side of the equation for Vj^q has the 
property that 

2it ^ 

U„ dT . 0, 


so that Vj^Q tends to zero exponentially when n The dependence of 

VjyjQ on the slow variable $ Induces the normal velocity component 

A 

U _ = AA' U From the equation of continuity we have 




n 

/ V. 

0 


MO 


dn , 


(5.5) 


which tends to a function of t when n Thus there exists a weak outer 
potential flow down by the mean velocity field in the Stokes layer. The outer 
potential flow decays algebraically in the normal direction and has a normal 
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velocity which matches with (5.5). However, the matching requires an 

“5 /8 

azimuthal velocity field of order g ' in the Stokes layer which tends to a 
given function of t and $ when n This velocity field is an 

'eigensolution' driven by the pressure field Pj^q which is then determined by 
matching with the outer potential flow. This outer flow is determined by 

first noting that when n can be written as 

OO 

= A($)A'($) 1' (5.6) 


where 


- 27T 

„-7/4 „ / r ; -inx , , 

= -2 Tq^ / J e dq dx , 


0 0 


MO 


n = ±1 , ±3 , ±5 , • • • , 


and 


g^ = 0, n = 0, i2 , ±4 , • • • . 


In order to determine the outer flow we write 


, . 1^1 6 ^'® , 


SO that the radial and azimuthal variations are now on the same length scale. 
The radial and azimuthal velocity components then expand as 


0/0 1/ “ S'*" 

U = g"^^^(vco) I ^ + •••} 


TT r" 2/8. . V? f V inx 1 

V = - g (vo)) ^ { I e 


(5.7) 
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where '1' =0 for n even whilst for n odd ¥ satisfies 

n n 


2 2 
S '? 3 "P 

n . n 


= 0 


3? 3$‘ 


'P •> 0, C ■> “, Y 0, |3> 1 -V CO 

n ’ ’ n ’ 


(5.8) 




The solution of (5.8) is 


f 

n 





A^(0)d9 
2 2 ' 


(5.9) 


The amplitude function A($) is determined as a solvability condition on the 

—3 /8 

differential system obtained by equating fundamental terms of order 3 
after substituting for U, V, W, and P from (5.2) into (2.5). The required 
condition is 


d^A r— »2-|, . 3 

— ^ + p[T - $ JA = yA 

do 


(5.10) 


where p is given by (2.15) with k = and ’’’O - '^Oc- The constant y is 
defined by 


T = 


2tt 00 

- / / {u P + V Q}dn dT 

0 0 


1 /'TT 00 J 

2^^* -ro'2 / / [v„ ■f- Idn dT 


3n 


V i 1 

C-^ 1' 


where 
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A A 


•|/\ AA AA A-IAA 

P = i k T- [V_ V_ + 2V„- V„1 - 3k U. U_ + U„ U_ 

2 c Oc'- 0 2 MO QJ c‘-2 0 2n 2 Ori- 


A A -1 A A A A I A A 

!u„ u_ + ± u_ U_ - U. U- - i. u_ u„ 

“ 2 Onpri 2 2n Onn On 2nn 2 0 2nnn- 


and 


A A 


A A 


A A 


Q == L"o ^2n + + 2V^ ^ VqJ 


The amplitude equation (5.10) must be solved subject to the condition 


A -> 0, 1$ 1 oo 

and of course reduces to (2.14) for A « 1. We postpone a discussion of the 

solution of (5.10) until after an investigation of the effect of a finite 

amplitude solution on the steady streaming of the basic flow. 

—3 /8 

The fundamental terms of order 6 ^ in (5.2) can be calculated when the 

solvability condition (5.10) is satisfied. The equations for the first and 

—3 /8 

second harmonic functions of order 3 can be solved directly without 

recourse to a solvability condition. The radial mean flow function Uj^q is 

—3 /8 

determined by (5.5) so that at order 3 it remains for us to discuss the 

azimuthal mean flow function . This function satisfies the equation 


9 _ V - 2 V = Am 
2 '^Ml 9 t Ml 


9n 


AA AA — '5//l ^0 

d« ^3, ("O h + “l ’o) * V ’ol • 


which is to be solved subject to 


''mi - "• 


q = 0,» 


(5.12) 
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However, the form of the nonlinear terms in (5.11) means that has a 
steady term in its Fourier series expansion. If the steady part of Is 
denoted by then the appropriate boundary conditions for are 



Thus the steady flow in the boundary layer induced by the finite amplitude 
disturbance does not decay to zero when n -^ «>. In fact we see from (5.11) 
and (5.12) that when p ->■ “ 


dA 


d$ 


where d is a constant to be calculated numerically. We found that 
d = -4.39 so that the azimuthal velocity component of the disturbance tends 
to -4.39 Uq P A when p ^ It is known (see Stuart [1966]) that 

the steady part of the azimuthal velocity component of the basic flow tends to 
3.2“^'^^ 3 - 3/8 ^ p ^ 00 with [6 - tt/2] = We see then 

that the steady streaming of the two-dimensional flow is modified by the 
instability. Moreover, it follows that in the outer steady streaming boundary 


layer the steady part of the basic flow and the instability cannot be found 

- 1 /, 

independently. This outer layer is of thickness aP ^ and if we take the 

variable K - (-^ ^)P we look for an outer steady flow given to first 

Si 


order by 






where 'P satisfies 
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\i/ =11/ vi; «. w w 


(5.14) 


which must be solved subject to 


->0, ? > 


'^Ocl-'^372 


1/2 ^ 

Oc d$J ’ 


(5.15) 


3$ _ 23/4 . /z A -^1 . C ^ 0. 


If we set A = 0 above we obtain the equations governing the attachment of a 

- 1/8 


steady streaming boundary layer within a 3 
In that case (5.14) is solved subject to 


neighborhood of 0 = it/2. 


= 0 , 


$ = 0 , 


(5.16) 


and the symmetry of (5.10) about 

dA 


applied since either A or 


do 


= 0 means that (5.16) can still be 

= 0 at 0=0. We note that for large 0 

the condition (5.15) reduces to ^ ~~3j^ , ^ ^ 0, so that, assuming that 

the boundary layer remains attached for finite values of 0 , the extra term 

dA 


proportional to A 


do 


in (5.15) merely produces an origin shift in the large 


dA 

the term A can be positive for some 


0 asymptotic solution of (5.14). However, we must recognize the fact that 

so that the slip velocity (5.15) 
can change sign at Intermediate values of 0 . If the magnitude of the 
inviscid slip velocity is sufficiently large where this occurs then the 
attached flow strategy falls and the steady streaming boundary layer will 
prematurely detach from the cylinder. This possibility does not occur for more 
general flows where the point of attachment of the steady streaming layer and 
the most unstable position do not coincide. In this case the steady streaming 
driven by the Instability is weak compared to that of the basic flow. 
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We return now to discuss the solution of (5.10) which of course depends 
crucially on the sign of j. Our calculations gave the totally unexpected 
result that 

Y = -.087, 

which means that finite amplitude solutions of (5.10) bifurcate subcritically 
from the eigenvalues of the linear problem and are unstable. In the work of 
Seminara and Hall [1976] the corresponding constant was found to be positive 
and this difference caused a great deal of concern. After exhaustive checks 
of the computer code, and in fact a repeat of the calculations using an 
independent program we believe that y is indeed negative and the 

bifurcations are Indeed subcritical. This could in principle be checked by a 
complete numerical investigation of the problem by Fourier expanding the 
velocity field in the z-direction but the complicated nature of the basic flow 
makes such a calculation non-trivial. 

6 . CONCLUSIONS 

We have shown that oscillatory viscous flows interacting with rigid 
boundaries of convex curvature can become unstable to Taylor-Gortler vortices. 
In particular, the flow Induced by the transverse oscillations of a circular 
cylinder is linearly unstable to Taylor-Gortler vortices localized where the 
slip velocity of the potential flow outside the boundary layer on the cylinder 
is a maximum. The results of our theory are in excellent agreement with 
Honji's observations over a wide range of values of the frequency parameter 
8 even though our results are formally valid only in the limit 3 

For an elliptical cylinder there are as yet no experimental results 
available. It would be interesting to see whether the cusp shaped curve for 
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a = 0 in Figure 4 could be found experimentally. There is no reason to 
suppose that the sensitive dependence of the critical Taylor number on the 
eccentricity and the angle of attack predicted in Section 4 could not he 
reproduced experimentally. 

The results of our nonlinear calculations are unexpected because of the 

prediction of the suhcritical nature of the instability. It is almost 

Invariably the case in the Taylor problem that the bifurcation to a Taylor 

vortex flow is supercritical but DlPrima and Sjbrand [1983] have found 

suhcritical bifurcation when considering the flow between counter-rotating 

cylinders. In fact, if we do not restrict the wavenumber to be that 

corresponding to the minimum on the neutral curve there will always be a 

finite band of wavenumbers where the Landau coefficient y is negative in the 

steady Taylor problem. This band of wavenumbers lies to the left of the point 

on the neutral curve where the wavenumbers on the left and right-hand branches 

are in the ratio 1:2. In the present problem the constant y becomes 

singular where the neutral values of the wavenumbers are K = .34 and 2K. 

In fact, near K the constant y behaves like - — so that to the left 

_ K - K 

of K there is a finite range of values of K for which y is positive. 

However, calculations show that the range of wavenumbers is only of length 

0(10“^) and y then becomes positive again. 

If the instability is indeed subcrltlcal then we presume that close to the 

critical Taylor number sufficiently large perturbations to the basic state 

will grow. It is possible that higher order nonlinear effects eventually 

cause these perturbations to equilibrate and that it is why Honji • observes 

some kind of steady state with Taylor-Gortler cells. In fact, even if 

nonlinear effects are not stabilizing at higher order then, because of the 

— 1 /fi 

localized nature of the Instability with the flow unstable in a 8 
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neighborhood of the most susceptible positions of the boundary layer we might 
expect that some periodicity along the cylinder would be observed. Indeed, it 
is known in parallel or nearly parallel flow stability theory that Tollmien- 
Schllchtlng waves can he observed even though they are subcritcally 
unstable. In the present problem, the subcritical nature of the bifurcation 
could be investigated by solving the full stability equations by Fourier 
analyzing in the z-direction and solving a large system of coupled nonlinear 
partial differential equations but such a computation would be nontrivial. 

Finally, we point out that perhaps Honji's results might in fact suggest 
that the instability does not develop supercritically in the manner usually 
found in the Taylor problem. We refer to the fact that Honjl gave two 
experimentally determined curves, one represents the onset of 'streaked flow' 
and a higher curve above which the streak could be observed because the flow 
was then separated and turbulent. We saw in the previous section that some 
finite amplitude solutions of (5.10) would cause the steady streaming boundary 
layer to separate prematurely. Thus our nonlinear calculations do in fact 
suggest an increasingly likely breakdown in the basic flow structure when the 
Taylor number is increased. Alternatively the separated flow observed by 
Honjl could be simply the unsteady two-dimensional separation of the Stokes 
layer on the cylinder. 

The author acknowledges some useful comments by Professor N. Riley in 
connection with the range of validity of the expansion procedures of Davidson 
and Riley for the elliptic cylinder problem. 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure Captions 


The neutral curve of the linear problem. 


The eigenfunctions U^, corresponding to the critical case 
Tq = 11.99, k = .51. 


A comparison between Honji's expermental points and linear theory. 


T 

The dependence of — on K for a = 0,0.2,0.4,it/2. 


The dependence of (j)^ on K for a = 0,0.2,0.4,tt/2. 


Figure 5. 
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